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Abstract. The third, fifth and sixth Painleve equations are studied by means of the 
weighted projective spaces CP 3 (p, q , r, s ) with suitable weights (p, q, r, s) determined by the 
Newton polyhedrons of the equations. Singular normal forms of the equations, symplectic 
atlases of the spaces of initial conditions, Riccati solutions and Boutroux’s coordinates are 
systematically studied in a unified way with the aid of the orbifold structure of C P 3 (p, q, r, s) 
and dynamical systems theory. 
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1 Introduction 

The first to sixth Painleve equations written in Hamiltonian forms are given by 

(P )• dx - dHj dy - dHj 

J ' dz dy ’ dz dx ’ 

J = I, II, IV, III(Z?8), III(Z?7), III(Z?6), V, VI, with the Hamiltonian functions defined as 

Hi = i x 2 - 2 y 3 - zy, 

Hu = i x 2 - ^ zy 2 - ay, 

Hjy = —xy 2 + x 2 y — 2 xyz — 2 ax + 2 j3y, 

TT 2 2 z I 

zH m( Ds ) = x y - — - -y, 

zH m(DT) = x 2 y 2 + zx + y + axy, 

zH m ( D6 ) = x 2 y 2 - xy 2 + zx + (a + f3)xy - ay, 

zHy = x(x + z)y(y — 1) + a 2 yz - a 3 xy - a\x(y - 1), 

z{z - l)Hyi = y(y -l){y - z)x 2 + a 2 {a\ + a 2 ){y - z) 

- (a 4 (y - 1 )(y - z) + a 3 y{y - z) + a 0 y{y - 1)) x, 

where ai,/3 G C are arbitrary parameters. Parameters of Hy i satisfy the constraint «o + a i + 
2a 2 + 0:3 + q ?4 = 0. The third Painleve equation is divided into three cases III(Zj>s), III(L>7), 
III(He) due to the geometry of the spaces of initial conditions [7, 8]. See [9] for the list of 
Hamiltonians and Backhand transformations written in these coordinates, although Hm^^ 
here is obtained by putting xm- l/(2y) for that in [9]. 

Among these Hamiltonians, only (Pi), (P11) and (Piv) are polynomials with respect to both 
of the dependent variables x, y and the independent variable 2. In [2, 3], (Pi), (P11) and (Piv) 
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are studied by means of a weighted projective space C P 3 (p,q,r,s), whose weight ( p,q,r,s ) is 
one of the invariants of the equation determined by the Newton polyhedron. In particular, the 
Painleve property, the spaces of initial conditions and Kovalevskaya exponents are investigated 
in detail. The purpose in this paper is to extend the previous result to the third, fifth and sixth 
Painleve equations, whose Newton polyhedrons are degenerate. 

According to [2, 3], let us recall the definition of the Newton diagram of a polynomial 
differential system. Consider the system of polynomial differential equations 


dxi 

dz 


fi(x i,.. .,x m ,z), 


i = 1,..., m. 


( 1 . 1 ) 


The exponent of a monomial x^ 1 included in the right hand side ft is defined as 

(pi,..., Hi-i, Pi — 1, fii+i, • • •, Umi d + 1). Each exponent specifies a point of the integer lattice 
in M m+1 . The Newton polyhedron of (1.1) is the convex hull of the union of the positive 
quadrants M+ +1 with vertices at the exponents of the monomials which appear in the system. 
The Newton diagram of the system is the union of the compact faces of its Newton polyhedron. 

We also consider the perturbative system 


= fi(x l, ■ •. ,x m ,z) + gi(x i, ...,x m ,z), i = 

where /,; and gi are polynomials. We suppose that 


( 1 . 2 ) 


(Al) the Newton polyhedron of the truncated system (1.1) has only one compact face and all 
exponents of monomials included in (1.1) lie on the face. 

In this case, there is a tuple of relatively prime positive integers (pi,... , p m ,r, s) and a hy¬ 
perplane p\Xi + • • • T PmX m + rz = s in M m+1 such that all exponents lie on the plane; i.e., any 
monomials x^ 1 ■ ■ ■ z v included in _/) satisfy 


Pi Pi H- Vpiim - 1) H- \-PmPm + r(r] + 1) = s. 


In [2], we further suppose that s — r = 1 though it is not so essential. To regard gt as a pertur¬ 
bation, we suppose that 

(A2) any monomials x^ 1 • • • x^z* 1 included in gt satisfy 

PiPi H - \~Pi(Pi — 1) H - 1 -PmPm. + r{g + l) < s 


(this implies that the exponents of gi lie on the lower side of the hyperplane). 

Due to the property of the Newton diagram, it is easy to verify that the truncated system (1.1) 
is invariant under the Z s -action given by 

(xi,.. • ,x m ,z) (w Pl xi,... ,u Prn x m ,uj r z), u:=e 2m/s . (1.3) 

We require the same for equation (1.2): 

(A3) Equation (1.2) is invariant under the Z s -action (1.3). 

A tuple of positive integers (p±,... ,p m ,r, s) is called the weight of the system (1.2). It is 
known that there is a one-to-one correspondence between nondegenerate Newton diagrams and 
toric varieties. If exponents he on the unique plane p\X\ + - ■ ■+p m x m +rx = s (assumption (Al)), 
then the associated toric variety is the weighted projective space CP m+1 (pi,... ,p m r, s), which 
is an (m + l)-dimensional orbifold, see Section 2.1 for the definition. 
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The first, second and fourth Painleve equations satisfy the above assumptions. For the first 
Painleve equation x' = 6y 2 + z, y' = x, put / = (6 y 2 + z, x ) and g = (0, 0). The Newton diagram 
is determined by three points (—1,2,1), (—1,0,2) and (1,-1,1). They lie on the unique plane 
3x + 2y + 4z = 5. For the second Painleve equation x' = 2 y 3 + yz + a, y' = x with a parameter a, 
put / = (2 y 3 + yz,x ) and g = (a,0). The Newton diagram is determined by points (—1,3,1), 
(—1,1, 2), (1, —1,1), which he on the plane 2x + y + 2z = 3. For the fourth Painleve equation, 
put / = (— x 2 + 2 xy + 2 xz, —y 2 + 2 xy — 2 yz) and g = (—2/1, —2a). The Newton diagram is given 
by the unique face on the plane x + y + z = 2 passing through the exponents ( 1 , 0 , 1 ), ( 0 , 1 , 1 ) and 
(0,0, 2). Hence, the weighted projective spaces associated with them are given by CP 3 (3, 2,4, 5), 
CP 3 (2, 1, 2,3) and CP 3 (1,1,1, 2), respectively. Note that g consists of terms including arbitrary 
parameters a, /3 for these systems. 

The orbifold CP m+1 (pi,... ,p m ,r, s) is regarded as a compactification of the phase space 
C m+1 = 

,x m ,z)} of the system (1.2). In [2, 3], the system (1.2), in particular the first, 
second and fourth Painleve equations, are studied with the aid of the geometry of CP m+1 (p \,..., 
p m ,r,s). In this paper, the third, fifth and sixth Painleve equations will be investigated, for 
which the Newton polyhedrons are degenerate and do not satisfy (Al). 

The third Painleve equation of type D§ is explicitly given by 


, . [zx' = -2x 2 y + 2 xy - (a + (3)x 4 

( III(D6)j: \ Z y> = 2xy 2 -y 2 + z + (a + P)y. 


(1.4) 


Put / = (—2 x 2 y + 2 xy, 2xy 2 — y 2 + z) and g = (—(a + f3)x + a, (a 4- /3)y). Exponents of / are 
given by (1,1, 0), (0,1,0), (0, —1,1). The Newton polyhedron generated by them has no compact 
faces because the positive quadrant with a vertex at (0,1,0) includes (1,1,0); the Newton 
diagram is empty. Nevertheless, these three points lie on the plane y + 2z = 1. Hence, we define 
the weight of (Pih(d 6 )) by ( 0 , 1 , 2 , 1 ) and use the weighted projective space CP 3 ( 0 , 1 , 2 , 1 ), which 
is not compact because of the nonpositive weight. 

The third Painleve equation of type Dj is given by 


(Piii(d 7 )) : 


zx' = —2 x 2 y — 1 — ax , 
zy' = 2 xy 2 + z + ay. 


Put / = (— 2x 2 y — 1, 2 xy 2 + z) and g = (—ax, ay). Exponents of / are given by (1, 1, 0), (—1, 0, 0) 
and (0, —1,1). The Newton diagram is again empty, however, these three points he on the plane 
—x + 2y + 3z = 1. Thus we define the weight of (Piii(d 7 )) by (—1, 2, 3,1) and use the weighted 
projective space CP 3 (—1, 2, 3,1). 

The third Painleve equation of type D$ is given by 


( P III(D 8 )) : 



1 9 

- 2 x 2 y 

2 y 

2 xy 2 . 


z 

V’ 


There are no parameters and put g = (0,0). Exponents of the system are (1,1,0), (—1,0,0) 
and (—1, —2,1). The Newton polyhedron is degenerate, however, these three points he on the 
plane —x + 2y + Az = 1. Thus we define the weight of (Pm(£> g )) to be (—1, 2,4,1) and use the 
weighted projective space CP 3 (— 1, 2,4,1). 

The fifth Painleve equation is given by 


(Pv): 


zx' = —2 x 2 y + x 2 + xz — 2 xyz + (a\ + a^)x — a^z, 
zy' = 2 xy 2 - 2xy - yz + y 2 z - (oq + a 3 )y + a\. 
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Put / = (—2 x 2 y + x 2 + xz — 2 xyz, 2 xy 2 — 2xy — yz + y 2 z) and g = ((oq + a^)x — z, —(a i + 
o.:i)y + a\). Exponents of / are (1,1,0), (1,0,0), (0, 0,1) and (0,1,1). Although there are four 
exponents, they lie on the unique plane x + z = 1. Thus we define the weight of (Py) by 
(1, 0,1,1) and use the weighted projective space CP 3 (1,0,1,1). 

The sixth Painleve equation is given by 

f z(z — l ) x ' = — 3x 2 y 2 + 2x 2 y + 2x 2 yz — x 2 z + g\, 

I z(z — 1 )y' = 2xy 3 — 2xy 2 — 2xy 2 z + 2 xyz + 52 , 

where ( 51 , 52 ) consists of terms including parameters. Exponents of the other terms are given 
by (1,2,0), (1,1,0), (1,1,1) and (1,0,1), which lie on the unique plane x = 1. Thus we define 
the weight of (Pyi) as (1, 0, 0,1) and use the weighted projective space CP 3 (1, 0, 0,1). 



(p, q, r , S ) 

deg(PI) 

K 

Ai, A2, A3 

l 

c 

Pi 

(3,2,4,5) 

6 

6 

6,4,5 

1 

2 

Pii 

(2,1,2, 3) 

4 

4 

4,2,3 

2 

3 

Piv 

( p l,1,2) 

3 

3 

3,1,2 

3 

4 

P III(D 8 ) 

(-1,2,4,1) 

2 

2 

2,4,1 

2 

3 

P III(D 7 ) 

(-1,2, 3,1) 

2 

2 

2,3,1 

2 

3 

P III(D 6 ) 

(0,1, 2,1) 

2 

2 

2,2,1 

3 

4 

p v 

(1,0,1,1) 

2 

2 

2,1,1 

4 

5 

p vi 

(1,0,0,1) 

2 

2 

2,0,1 

5 

6 


Table 1. deg(R) denotes the weighted degree of the Hamiltonian function with respect to the weight 
deg(cc) = p, deg (y) = q , deg(^) = r. n denotes the Kovalevskaya exponent defined in Section 2.2. 
(Ai, A 2 , A 3 ) is the weight for the weighted blow-up, which coincides with the eigenvalues of the Jacobi 
matrix at the singularity. In [2], it is proved that deg (H) = k = Ai and r = A 2 , s = A 3 . I gives the 
number of types of Laurent series solutions given in Section 2.2, and c is the number of local charts for 
the space of initial conditions. We always have c = l + 1. 

In this paper, the third, fifth and sixth Painleve equations are studied with the aid of the 
weighted projective spaces CP 3 (p, 5 , r, s ) and dynamical systems theory. The weighted projective 
space is decomposed into the disjoint sum CP 3 (p, q, r, s) ~ CP 2 (p, q, r) U C 3 . This implies that 
the natural phase space C 3 = {(&, y, z)} of the Painleve equations is embedded in CP 3 (p, q , r, s ) 
and the two-dimensional manifold C P 2 (p,q 1 r) is attached at infinity. We regard the Painleve 
equation as a three-dimensional autonomous vector field defined on C P 3 (p, q,r, s). Then, the 
vector field has several fixed points on the infinity set C P 2 (p,q,r). These fixed points de¬ 
scribe the asymptotic behavior of solutions. Some of these fixed points correspond to movable 
singularities, and the other correspond to the irregular singular point. 

The dynamical systems theory is applied to these fixed points to investigate the Painleve 
equations. The singular normal form of the Painleve equation [3, 6 ], which is a local integrable 
system around a movable singularity, is obtained by applying the normal form theory around 
fixed points. The space of initial conditions and its symplectic atlas are constructed by the 
weighted blow-up at these fixed points. The weight for the weighted blow-up, which is also an 
invariant of the Painleve equation related to the Kovalevskaya exponent [2, 3], is determined by 
eigenvalues of the Jacobi matrix of the vector field at the fixed points. It is known that the 
Painleve equations are reduced to the Riccati equations when the parameters take certain specific 
values. Such Riccati solutions are characterized as a center (un)stable manifold at the fixed point 
on CP 3 (p, q , r, s ). Although some of these results are well known for experts, our new approach 
based on the weighted projective space and dynamical systems theory provides a systematic way 
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to investigate them. From our analysis, it turns out that the weights and the Kovalevskaya expo¬ 
nents are important invariants of the Painleve equations. In particular, the Painleve equations 
may be classified by these invariants, which will be reported in a forthcoming paper. 

Our method will be explained in detail for the third Painleve equation of type D§ in Section 3. 
Since the strategy for the other Painleve equations (Pmm 7 )), (Piii(d 8 )), (Py) and (Pyi) is 
completely the same as that for (Phi(d 6 )), we only show a sketch and several formulae for them 
after Section 4. See [3] for (Pi), (Pn) and (Piv). 

2 Settings 

2.1 Weighted projective spaces 

For a tuple of integers (p\,... ,p m , r, s ), consider the weighted C*-action on C m+2 defined by 
(x \,.. .,x m ,z,£) (A Pl zi,..., X Pm x m , X r z, X s e), A e C* := C\{0}. 

The quotient space 

C P m+1 (pi, ...,p m ,r, s ) := C m+ 2 \{ 0 }/C* 

gives an (m + l)-dimensional orbifold called the weighted projective space. In this paper, we 
only use a three-dimensional space. When ( p , q , r, s) are positive integers, the orbifold structure 
of C P 3 (p, q , r, s) is obtained as follows: 

The space CP 3 (p, q, r, s ) is defined by the equivalence relation on C 4 \{0} 

(x, y, z, e) ~ ( X p x , X q y, A r z, X s e). 

(i) When x ^ 0, 

(x,y,z,e) ~ (1 ,x~ q/p y,x~ r/p z,x~ s/p £) =: (1, Y u Z u ei). 

Due to the choice of the branch of x 1 / p , we also obtain 

(Ki,Zi,ei) ~ (e- 2 ^Ij,e- 2 ™/^i,e- 2 s "/p £ i), 

by putting x i —> e 2m x. This implies that the subset of C P 3 (p,q,r,s) such that x ^ 0 is homeo- 
morphic to C 3 /Z p , where the Z p -action is defined as above. 

(ii) When y ^ 0, 

(x, y, z, e) ~ ( y~ p/q x , 1, y~ r/q z, y~ s,q e) =: (X 2 , 1, Z 2 ,e 2 ). 

Because of the choice of the branch of y 1 ^ q , we obtain 
(X 2 ,Z 2 ,s 2 ) ~ (e- 2p ™/ q X 2 ,e- 2r ™/ q Z 2 ,e- 2s7Ti / q e 2 ). 

Hence, the subset of C P 3 (p, q, r, s ) with y ^ 0 is homeomorphic to C 3 /Z g . 

(iii) When z/0, 

{x,y,z,e) ~ (z~ p/r x,z~ q/r y, l,z~ s/r e) =: (X 3 , Y 3 , 1, e 3 ). 

Similarly, the subset {z / 0} C C P 3 (p, q, r, s) is homeomorphic to C 3 /Z r . 

(iv) When e / 0, 

(x,y,z,s) ~ (e~ p/s x,£~ q/s y,£~ r/s z,l) =: (X 4 , Y 4 , Z 4 , 1). 

The subset {e / 0} C C P 3 (p, q, r, s) is homeomorphic to C 3 /Z s . 
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This proves that the orbifold structure of C P 3 (p,q,r, s ) is given by 

C P 3 (p, q, r, s ) = C 3 /Z p U C 3 /Z q U C 3 /Z r U C 3 /Z s . 

The local charts (Yi,Zi,£i), (A' 2 , Z 2 , e 2 ), (^ 3 ,^ 3 , £ 3 ) and {X^Y^ Z 4 ) defined above are called 
inhomogeneous coordinates as the usual projective space. Note that they give coordinates on 
the lift C 3 , not on the quotient C 3 /Zj (i = p,q,r,s ). Therefore, any equations written in these 
inhomogeneous coordinates should be invariant under the corresponding Zj actions. 

In what follows, we use the notation (x, y, z) for the fourth local chart instead of (X 4 , Y 4 , Z 4 ) 
because the Painleve equation will be given on this chart. The transformations between inho¬ 
mogeneous coordinates are give by 

x = = X*? 1 ' = X 3 ej r/ ‘, y = Y,Ep'- = Ep'- = Y 3 e 3 ,, ‘, 

Z = Z 1 E, r/, = Z 2 E? /S =E^ /S . (2.1) 

The same transformation rule holds even if p , q, r, s include negative integers. If there are 0 
among them, for example if p = 0, then we have CP 3 (0, q,r,s) — C x C P 2 (q, r, s). CP 3 (p, q, r, s ) 
is compact if and only if all p, q, r, s are positive. 

2.2 Laurent series solutions and Kovalevskaya exponents 

To construct the space of initial conditions, we need the expressions of the Laurent series of 
solutions. Let (p, q , r, s ) be the weight of a given system determined by the Newton polyhedron. 
Suppose that the system has a Laurent series solution of the form 

OO OO 

x = ^ - zo)~P +n , y = J 2 B n {z - zo)~ q+n , 

n =0 n =0 

where (^4o>-E>o) / (0,0) and zq is a movable pole. Such a Laurent series solution is called 
regular. A Laurent series solution is called exceptional if it is not expressed in this form; i.e., 
(Aq,Bq) = (0,0) or the order of a pole of either x or y is larger than p or q, respectively. If 
a regular Laurent series represents a general solution of the system, it includes an arbitrary 
parameter, which depends on initial conditions, other than zq. The smallest integer n such that 
(A k , B k ) includes an arbitrary parameter is called the Kovalevskaya exponent. In [2], it is proved 
that the Kovalevskaya exponent of the regular Laurent series solution is invariant under a certain 
class of coordinates transformations including the automorphism group of C P 3 (p, q,r, s). For 
the first, second and fourth Painleve equations, all Laurent series solutions are regular because 
they satisfy the assumptions (Al) and (A2). The third, fifth and sixth Painleve equations have 
exceptional Laurent series solutions, however, they can be converted into the regular series by 
the Backlund transformations. Hence, the Kovalevskaya exponents k of exceptional Laurent 
series solutions are well-defined and given as in Table 1. In what follows, denote T := z — zq. 

(Piii^e)): The third Painleve equation of type D§ has three types of Laurent series solutions 
given by 

(i) x = 1 + ^--T + A 2 -T 2 + 0(T 3 ), 

zo 

y = -z 0 • T ” 1 + ^(-1 - a + /3) + B 2 ■ T + 0(T 2 ), 

(ii) x = 0-T° - ^-T + A 2 -T 2 + 0(T 3 ), 

Zo 

y = z 0 • T~ l + i(l - a + $) + B 2 ■ T + 0(T 2 ), 
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(hi) x = -z 0 -T ~ 2 + 0-T - 1 +A 2 + 0(T), 

y = —T + ~ 2 + a + f 3 . r 2 + S 2 • T 3 + 0 (T 4 ), 
izo 

where A 2 is an arbitrary constant and B 2 is a certain function of A 2 . Since (p,q) = (0,1), the 
first two series are regular, while the last one is exceptional. The Kovalevskaya exponents of all 
series are k = 2. 

(Piii(d 7 )) : The third Painleve equation of type D 7 has two types of Laurent series solutions 
given by 

(i) x = -■T + ^f^-T 2 + A 2 -T 3 + 0{T 4 ), 

z 0 ^z 0 

y = -4 • T' 2 - z o • T~ l + B 2 + O(T), 

(ii) x = —zo ■ T ~ 2 + 0 • T ~ 1 + A 2 + 0(T), 

y = —T 1 + • T 2 + B 2 ■ T 3 + 0(T 4 ), 

2 zo 

where A 2 is an arbitrary constant and B 2 is a certain function of A 2 . Since (p. q) = (—1,2), the 
former series is regular, while the latter one is exceptional. The Kovalevskaya exponents of both 
series are k = 2. 

(Piii(d 8 )) : The third Painleve equation of type Dg has two types of Laurent series solutions 
given by 

(i) x = --j- • t + -L • r 2 + 2B * ~ 1 • r 3 + o(t 4 ), 

2 Z 0 4 zl 4 zq v ’ 

y = 2 zl • T- 2 + 2zo ■ T _1 + B 2 + 0(T), 

(ii) x = 2 zl- T ~ 3 + 2z 0 • T” 2 + 0 • T~ l + 0(1), 

y = —T 2 + 0 • T 3 + B 2 ■ T 4 + 0(T 5 ), 

2^o 

where B 2 is an arbitrary constant. Since (p,q) = (—1,2), the former series is regular, while the 
latter one is exceptional. The Kovalevskaya exponents of both series are k = 2. 

(Py): The fifth Painleve equation has four types of Laurent series solutions given by 

(i) x = zq ■ T 1 +-(1 — 20 + ai — 03 ) + ^2 • T + O(T^) 1 

y = l + — - T + B 2 -T 2 + 0(T 3 ), 
zo 

(ii) x = — zo ■ T 1 + - (—1 — zo + aq — 0 : 3 ) + A 2 ■ T + O (T 2 ), 

y = O.T 0 - — - T + B 2 - T 2 + 0(T 3 ), 

zo 

(hi) x = —zo + {011 + o 2 + 03 — 2 ) • T + A 2 • + 0(T 4 ), 

V = T- ' + zo + + y* 2 + 03 ~ 2 + b 2 ■ t + o (r 2 ), 

2^0 

(iv) x = 0-T° + a 2 -T + 2 l 2 -T 2 + O(T 3 ), 

» = -T-' + 20 + ai + 202 + 03 + Eh ■ T + o(r 2 ). 


where A 2 is an arbitrary constant and 1? 2 is a certain function of 7l 2 . Since (p. q) = (1,0), 
(i) and (ii) are regular, while (iii) and (iv) are exceptional. The Kovalevskaya exponents of all 
series are k = 2. 



H. Chiba 


(Pvi): The sixth Painleve equation has five types of Laurent series solutions given by 


(i) 


(ii) 


(hi) 


(iv) 

(v) 


x = T _1 + 


2 — 4^0 + «0 — 2zo«0 + ^0^3 — 04 + ZQOL4 


2 z 0 {z 0 - 1 ) 
y = z o + (2 + ao) • T + B 2 • T 2 + 0(r 3 ), 


+ ^ 2 -t + o(t 2 ), 


x = -2 0 • T 1 + 


1 — zq — oiq + 2 a 3 — zqcx^ — <y,\ + zoQq 
2(2:0 - 1) 


+ a 2 -r + o(T 2 ), 


y = l-^-.T + B 2 -T 2 + 0(T 3 ), 

20 

/ rri —1 , z 0 - a 0 — z 0 a 3 + a 4 + Z0CX4 , 2 \ 

x = (z 0 -l)-T H---h • T + Oir z ), 

2zo 

y = 0 • T° + ° 4 •T + B 2 -T 2 + 0 (T 3 ), 

20-I 

ai(ai + a 2 ) , nt(rr 2 \ z 0 (z 0 - 1) x 

x =- - -— ■T + OyT), y= - T +0(1), 

2:0(20 - 1 ) ai 

CVl(% 2 rp r 4 ( r p‘2\ 2o(20 1) rri—1 , .^v /-i \ 

x = —, -—-r + o(T), y = - T +°i> 

20(20 - 1 ) ai 


where A 2 is an arbitrary constant and B 2 is a certain function of A 2 . Since (p,q) = (1,0), 
(i), (ii) and (iii) are regular, while (iv) and (v) are exceptional. The Kovalevskaya exponents of 
all series are k = 2 . 

For all Painleve equations, the number of types of Laurent series solutions is smaller than 
the number of local charts of the space of initial conditions by one, see Table 1. 


3 The third Painleve equation of type D 6 

3.1 Pni(z? 6 ) on CP 3 (0,1, 2,1) 

The orbifold structure of CP 3 (0,1, 2,1) is given by 

CP 3 ( 0 , 1 , 2 , 1 ) = C x CP 2 ( 1 , 2 , 1 ) = C x (c 2 u c 2 /z 2 u c 2 ) = c 3 u (c x c 2 /z 2 ) u c 3 . 

Thus, the space is covered by three inhomogeneous coordinates (X 2 , Z 2 , £2), (X 3 , Y3, £3) and 
( x , y, z) related as 

x = X 2 = X 3 , y = £2 1 = l3£ 3 \ 2 = Z 2 £2 2 = +C 2 (3-1) 

(the first local chart (Yf, Z\,£i) does not appear because p = 0). 

We give the third Painleve equation of type D$ on the local chart (x, y, z). On the other local 
charts, (Pni(n 6 )) is expressed as 

dX 2 = 2X2 - 2X 2 + £ 2 (q - aX 2 - /3X 2 ) 
de 2 £2(1 — 2X2 — Z 2 — ae 2 — (3e 2 ) 

dZ 2 _ 2Z 2 — 4X2^2 — 2Z$ + £2^2(1 — 2a — 2/3) 
de 2 £2(1 — 2X2 — Z 2 — ae 2 — /3e 2 ) 

^ = -4 (2X 3 K 3 - 2X|l3 + £ 3(0 - 0 X 3 - /3X 3 )) , 
as 3 s ^ 

S = “ lj i 1 - ** + 2X3132 + ^ (- \ + “ + P) ) ' 
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In order to apply dynamical systems theory later, it is convenient to rewrite them as 3-dimen¬ 
sional autonomous vector fields of the form 

X 2 = 2X 2 - 2Xl + £ 2 (« - al 2 - f3X 2 ), 

Z 2 = 2 Z 2 — AX 2 Z 2 — 2 Z\ + £ 2 ^ 2(1 — 2 a — 2/3), 

£2 = £ 2(1 — 2 X 2 — Z 2 — ae 2 — @£ 2 ), (3-2) 

and 

X 3 = 2 X 3 F 3 - 2 X%Y 3 + £ 3 (a - aX 3 - 0X 3 ), 

*3 = 1- *3 + 2 X 3 P 3 2 + £ 3 Y 3 (-l + a + /3), 

£3 = -f, (3.3) 

where (") = d/dt and t £ C parameterizes each integral curve. 

We also have the decomposition 

C-P 3 (0,1, 2,1) = CP 2 (0,1, 2) U C 3 (disjoint). 

Since CP 2 (0, 1 , 2 ) = C x CP 1 (1, 2 ) and CP 1 (1, 2 ) is the usual projective line, 

CP 3 (0,1, 2,1) = (C x CP 1 ) U C 3 , 

where C 3 = {(x, y, z )} and C x CP 1 = {{X 2 , Z 2 ,0)} U {(X3, 13 ,0)} in coordinates. This implies 
that the set C X CP 1 is attached at infinity of the natural phase space C 3 = {(x,y,z)} of 
(Phi(£) 6 )); and the asymptotic behavior of solutions can be studied by the limit £2 —> 0 or 
£3 ->• 0 . 

The vector fields (3.2) and (3.3) have exactly four fixed points on the infinity set {£2 = 
0} U {£3 = 0} given by 

(X 2 ,Z 2 ,£ 2 ) = (0,0,0), (1,0,0), (0,1,0), (1,-1,0). 

The Jacobi matrices of (3.2) at the fixed points are 

/2 0 cA (2 0 f3\ / 2 0 a\ / 2 0 f3\ 

0 2 0, - 0 2 0, -4 -2 7 , - -4 -2 7 , (3.4) 

\0 0 1 / \0 0 1 / \° 0 0 / \° 0 0 / 

respectively, where 7 = 1 — 2a — 2/3. The latter two points, for which Z 2 7 ^ 0 =$■ z = 00, 
correspond to the irregular singular point. Since the Jacobi matrix has a zero eigenvalue, there 
exists a one-dimensional center manifold at each fixed point. The asymptotic expansion of the 
center manifold [5] yields the asymptotic expansion of (x(z),y(z)) as \z\ —> 00. 

On the other hand, the former two fixed points correspond to movable poles. Indeed, it 
is easy to verify by using (3.1) that the regular Laurent series solutions (ii) and (i) given in 
Section 2.2 converge to the points (X 2 , Z 2 ,e 2 ) = (0,0,0) and (1,0,0), respectively, as z —> zq. 
The exceptional Laurent series solution (iii) does not converge to some point on CP 3 (0,1,2,1) 
as z —> zo ; the space CP 3 (0,1, 2 , 1 ) is not compact. 

To treat the Laurent series (iii), we use the Backlund transformations. See [9] for the complete 
list of the Backlund transformations of two-dimensional Painleve equations. It is known that 
the transformation groups of the Painleve equations are isomorphic to the extended affine Weyl 
groups. For a classical root system R , the affine Weyl group and the extended affine Weyl 
group are denoted by VF(i?^)) and W(R^), respectively. Let G = Aut(J?^^) be the Dynkin 
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automorphism group of the extended Dynkin diagram. We have W(R^) = G ix W(R (1) ). For 
the third Painleve equation of type Dq, R = 2A\ and 

W((2A 1 )«) = Gk W((2Ai)«), W((2A 1 )W) = < So ,si,s' 0 ,4), 

G = Aut (D^ ] ) = Aut (( 2 Ai) (1) ) = < 7 ri, 7t 2 , erf). 

The action of each element is given in Table 2, where /o = x + (a + (3 — 1 )/y + z/y 2 , f\ = 
Py + (x — 1 )y 2 + z, f 2 = ay + xy 2 + 2 and / 3 = 2 xy + a + P - 1. 



a 

P 

y 

X 

z 

So 

2 — a 

P 

1 — a 

(1 — a ){fo — 2 y) (1 -a) 2 z 

z 

V + fo ~ 1 

x h f 1 2 

Si 

—a 

P 

a 

y + - 

X 

X 

z 

4 

a 

2-P 

1 ~P 
y+ , 

Jo 

(1 — P)h (1 -P) 2 z 

X h /| 

z 

S| 

a 

~P 

«+ \ 

X 

z 

7Tl 

1 — a 

P 

z 

y 

^(xy-y + P) + 1 

z 

7T 2 

a 

1 ~P 

z 

y 

-^(xy + a) 

z 

°1 

P 

a 

-y 

1 — X 

—z 


Table 2. The action of the extended affine Weyl group for (Pih(_d 6 ))- 


Let us consider another space CP 3 (0,1,2,1) with inhomogeneous coordinates (X 2 , Z 2 ,e 2 ), 
(A 3 , Y 3 , e 3 ) and (x, y, z) satisfying the same formula as (3.1). We glue two copies of CP 3 (0,1,2,1) 
by the transformation n 2 


{x, y, z, a, (3) = n 2 (x, y, z, a, /?) = 


-~(xy + a),-,z,a, 1 -/3 ) , 

z y 


(3.5) 


which defines a manifold denoted by Adi. Since (Piii(_d 6 )) is invariant under the action of tt 2 , the 
system written in the {X 2 , Z 2 , ? 2 )-chart also satisfies equation (3.2), in which (a,/3) is replaced 
by (5, j3) = (a, 1 — /3). By the formula (3.5) or 

AI 2 = ~{xy 2 + ay)/z, Z 2 = y 2 /z, e 2 = y/z, 


it is easy to show that the exceptional Laurent series (iii) is converted to the regular series (i) 
in the (x, y, 2f)-chart, and it approaches to the fixed point (A 2 , Z 2 , e 2 ) = (1, 0, 0) as z —> zq. 

We have proved that Laurent series solutions (i), (ii) and (iii) of (Piii(_d 6 )) converge to the 
fixed points (X 2 , Z 2 , s 2 ) = (1, 0, 0), (0,0,0) and (X 2 , Z 2 ,s 2 ) = (1,0,0) on Adi, respectively, as 
z —>■ z o- Hence, the study of movable poles are reduced to the study of the fixed points and 
dynamical systems theory is applicable to investigate them. 

As a simple application, we can prove the next theorem. 


Theorem 3.1. There exists a local analytic transformation (X 2 , Z 2 ,e 2 ) (->■ ( u,v,w) defined near 
(X 2 , Z 2 , e 2 ) = (0,0,0) such that equation (3.2) is transformed into the linearized system 


u = 2 u + aw, v = 2 v, vj = w. 


Similarly, the vector field is locally linearized around (X 2 , Z 2 ,e 2 ) = (1,0,0) and {X 2l Z 2 ,e 2 ) = 

( 1 , 0 , 0 ). 
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This result implies that (Piii(d 6 )) is locally transformed to the integrable system around 
each movable singularity. A proof is a straightforward application of Poincare’s linearization 
theorem [5] in dynamical systems theory. See [3] for the detail, in which a similar result is 
proved for the first, second and fourth Painleve equations, and also [2], in which it is proved 
that any differential equations having the Painleve property is locally linearizable. Such a local 
integrable system for the Painleve equation is called the singular normal form in [3, 6 ]. 

3.2 The space of initial conditions 

A purpose in this section is to construct the space of initial conditions for (Pih(_d 6 )). On the 
manifold Xi\, there are three singularities of the foliation of integral curves; {X 2 ,Z 2 , £ 2 ) = 
(1, 0,0), (0,0,0) and (^ 2 ,^ 2 , £ 2 ) = (1,0,0). They correspond to movable poles of the Laurent 
series solutions (i), (ii) and (iii), respectively. We will resolve these singularities by weighted 
blow-ups. On the blow-up space, (Piii(d 6 )) is again written in a Hamiltonian system, whose 
Hamiltonian function is polynomial in dependent variables. This implies that singularities of 
the foliation are resolved and the space of initial conditions is obtained by three times blow-ups 
of All. This strategy is applicable to the other Painleve equations, even for higher-dimensional 
Painleve equations. 

(i) blow-up at (X 2 ,Z 2 ,S 2 ) = ( 0 , 0 , 0 ). 

For the vector field (3.2), put 

U = X 2 + CXS2 , v = Z 2 , w = £ 2 - 

Then, the linear part is diagonalized and we obtain 

u = 2u — 2u 2 + w(au — av — f3u), 
v = 2 v + v(w — Au — 2 v + 2 aw — 2 f3w), 
w = w — w(2u + v — aw + (3w). 

The origin ( u , v, w ) = (0, 0, 0) is a singularity of the foliation of integral curves. To resolve it, 
we introduce the weighted blow-up defined by 

U = lt\ = v\u 2 = W^lls, V = u\v 1 = v\ = w\v 3, W = U\W\ = V2W2 = W3. 

The weight (2, 2,1), the exponents in the right hand sides, is taken from the eigenvalues of the 
Jacobi matrix at the singularity. The exceptional divisor {u\ = 0} U {v 2 = 0} U {u >3 = 0} is the 
weighted projective space CP 2 (2,2,1). The relation between the original coordinates ( x,y,z ) 
and the new coordinates (us,vs,ws) is given by 

x = u 3 wl - aw 3 , y = l/w 3 , z = v 3 , 

u 3 = xy 2 + ay, w 3 = 1/y, v 3 = z. (3.6) 

Note that eventually the independent variable 2 is not changed and (3.6) defines a fiber bundle 
over z-space, whose fiber is a (x, y)- space and (u 3 ,w 3 )-sp&ce glued by the above relation. In the 
new coordinates, (Piii(d 6 )) is transformed into the Hamiltonian system 

dw 3 dH\ du 3 dH\ 

dz du 3 ’ dz dw 3 ’ 

with the Hamiltonian function 


zH\ = —u + u 2 w 2 + uw 2 z — awz — (a — /3)uw 
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(here, the subscript is omitted for simplicity). Since H\ is polynomial in (113,^3), there are 
no singularities of the foliation in this chart; the singularity associated with the Laurent series 
solution (ii) is resolved. Furthermore, we can verify that 

dx A dy = dw 3 A du 3, dx A dy + dH A dz = dw 3 A d«3 + dffi A dz. (3.7) 

(ii) blow-up at ( X 2 ,Z 2 ,£ 2 ) = (1,0,0). 

Putting X 2 — 1 = A2 and u = A^2 + @£ 2 , v = Z2, w = £2 for the vector field (3.2) results in 

h = —2u — 2 u 2 — u;(cra — (3u + /3u), 
u = — 2v + v(w — 4 u — 2u — 2cra + 2 /din), 
re = —re — w{2u + v + aw — /3rc). 

To resolve the singularity at ( u,v,w ) = (0,0,0), we introduce the weighted blow-up with the 
weight (2, 2,1) 

U = u\ = U5U5 = WqU6, V = u\v 4 = V§ = WqV 6, U? = «4U;4 = V5W5 = Wq. 

The relation between the original coordinates ( x,y,z ) and the new coordinates (uq,vq,wq) is 
given by 

x = uqwI - / 3 we + 1, y = l/we, z = v&, 

u 6 = xy 2 + [ 3 y - y 2 , w 6 = l/y, v 6 = z, (3.8) 

which defines a fiber bundle over 2-space. In the new coordinates, (Pih(d 6 )) written as the 
Hamiltonian system 

dw& dH2 due dH2 

dz due ’ dz Owe ’ 

with the Hamiltonian function 

ZH2 = u + u 2 w 2 + uw 2 z — ( 3 wz + (a — / 3 )uw (3-9) 

(here, the subscript is omitted for simplicity). Since H2 is polynomial in (ue, we), the singularity 
associated with the Laurent series solution (i) is resolved. As before, we can verify the symplectic 
relation (3.7). 

(iii) blow-up at (X 2 , £ 2 ,^ 2 ) = (1,0,0). 

This singularity is resolved by the same way as above by the weighted blow-up with the 
weight (2,2,1). The final result is easily obtained by the Backlund transformation (3.5) as 
follows. Define the new coordinates (ug,vg,wg) by 

ug = xy 2 + fiy — y 2 , w 9 = 1/y, v 9 = z, 

which is the same formula as (3.8). Substituting (3.5) yields 

ug = -z 2 /y 2 - xz - (a + /3 - l)z/y, w 9 = y/z, vg = z. (3.10) 

By this coordinates change, (Piii(_d 6 )) is transformed into the Hamiltonian system of (ug,wg), 
whose Hamiltonian function has the same form as (3.9), though (a,/3) is replaced by (a, /3) = 
(a, 1 - / 3 ). 

In this manner, all singularities are resolved and we have 

Theorem 3 . 2 . The space of initial conditions E(z) of { Pih(_d 6 )) given by C 2 xy ^ u c («,«> u 
C(u 6 «, 6 ) u C ( U9 u , (J ) glued by the symplectic transformations (3.6), (3.8) and (3.10). The space E{z) 
is a nonsingular symplectic surface parameterized by z G C\{0}, on which (Phi(_d 6 )) is expressed 
as a polynomial Hamiltonian system. 
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3.3 The Riccati solutions 

It is known that when the parameter a is an integer or ft is an integer, there exists a one- 
parameter family of solutions of (Pmm 6 )) satisfying the Riccati equation, which is equivalent to 
the Bessel equation. For example, when a = 0 (resp. j3 = 0), the Riccati equation is given by 

dy 2,o dy 2 i i /q 11 \ 

z—~ = —y + py + z, resp. z— = y + ay + z. (3.11) 

az az 

The Backlund transformations yield the Riccati equations for a general case a 6 Z or ft E Z. 
Let us prove this fact from a view point of dynamical systems theory. 

Recall that there are two fixed points (X 2 , ^ 2 ,^ 2 ) = (0,1,0), (1, —1, 0) of the vector field (3.2) 
corresponding to the irregular singular point. The Jacobi matrices at these points are shown 
in (3.4), and eigenvalues of them are 2, —2, 0. Hence, there exist a center-stable manifold and 
a center-unstable manifold at these points. Let us calculate them explicitly. 

For the point (X 2 , Z 2, £2) = (0,1, 0), put Z 2 = Z 2 — 1 and 

u = X 2 +Z 2 +^£ 2 ( 01 + 2/3 - 1), v = —X 2 - ^a£ 2 , w = £ 2 . 

Then, the linear part of equation (3.2) is diagonalized around (0,1,0) and we obtain 


u = —2 u + u ^—2 u — -w — — w + fiiv'j + (1 — a) + - w 2 — 

v = 2v + v ^2v + —w — /3w^j + -a ^ uw + -w 2 — (3w 2 ^j , 


w = w | v — u — -w 


The stable, unstable, center subspaces are u, v, re-directions, respectively. In particular, the 
center-unstable manifold (resp. the center-stable manifold) is expressed as the graph of some 
function v = 4 >(u,w) (resp. u = ip(v,w)). 

It is obvious that when a = 0, the center-unstable manifold is exactly given by v = 0. 
This implies X 2 = x = 0. When a = 0 and x = 0, equation (1.4) is reduced to the Riccati 
equation (3.11). Similarly, when a = 1, the center-stable manifold is exactly given by u = 0. 
This implies 


0 — X 2 + Z 2 + @£2 — x + z/y 2 + P/y — 1 . 

Substituting this relation to equation (1.4) yields the Riccati equation for a = 1, though it is 
also obtained by the Backlund transformation to that for a = 0. The same argument at the 
point (X 2 , Z 2 ,£ 2 ) = (1, —1, 0) provides the Riccati equation for fi = 0. 

Theorem 3.3. When a G Z or (3 G Z, there exists a one-parameter family of solutions of 
(Piii(d 6 )) governed, by the Riccati equation of Bessel type. The one-parameter family of solutions 
forms a center-(un)stable manifold of (3.2) in CP 3 (0,1,2,1). 

The same argument is applicable to the other Painleve equations, even for higher-dimensional 
Painleve equations. For two-dimensional equations, it is well known that when parameters take 
certain specific values, the Painleve equations are reduced to Riccati-type equations except 
for (Pi), (Piii(d 7 )) an d (Piii(d 8 ))- F° r such specific values of parameters, center-(un)stable 
manifolds in C P 3 (p, q,r, s) are exactly calculated and the Riccati equations are obtained by 
restricting the equations to the center-(un)stable manifolds. 



14 


H. Chiba 


The result is summarized in Table 3. The third column shows one of the parameters for which 
equations are reduced to the Riccati equations. The other possible parameters are obtained by 
the Backhand transformations. The fourth column denotes the name of the Riccati equation 
when it is written as a second order linear equation. Note that the weight of the Riccati equa¬ 
tion is also defined through the Newton polyhedron. For example, the Riccati equation of Airy 
type is defined by 


dy 

dz 


y 2 + z. 


The exponents of monomials in the right hand side are (1,1) and (—1, 2 ). Since they are on the 
line y+2z = 3, the weighted projective space for the equation is CP 2 ( 1 , 2,3). See [ 4 ] for the analy¬ 
sis of the Airy equation by means of the weighted projective space. The last column of Table 3 
gives the weight of each Riccati equation. It is interesting to note that these weights are obtained 
by deleting the first or second numbers from the weights of the corresponding Painleve equations. 



weight 

parameter 

Riccati 

weight 

Pll 

(2,1,2, 3) 

a = 1/2 

Airy 

(1,2,3) 

Piv 

( 1 , 1 , 1 , 2 ) 

a = 0 

Hermite 

( 1 , 1 , 2 ) 

PlIICDe) 

( 0 , 1 , 2 , 1 ) 

a = 0 

Bessel 

( 1 , 2 , 1 ) 

Pv 

( 1 , 0 , 1 , 1 ) 

ai = 0 

CHG 

(M,l) 

Pvi 

( 1 , 0 , 0 , 1 ) 

ai = 0 

HG 

( 1 , 0 , 1 ) 


Table 3 . The type of Riccati equations and their weights. CHG and HG denote the confluent hyper¬ 
geometric and hypergeometric equations, respectively. 


3.4 Boutroux’s coordinates 

For the first and second Painleve equations, the third local chart (A 3 , 5^3, £3) of C P 3 (p,q,r,s) 
defined by equation (2.1) is equivalent to Boutroux’s coordinates introduced in [1] to investigate 
the irregular singular point 2: = 00 . For the other Painleve equations, the chart (A3, 13 , £3) plays 
the same role as Boutroux’s coordinates; as £3 — > 0, equation (3.3) is reduced to the autonomous 
Hamiltonian system 

A 3 = 2A3I3 - 2 Af Y3, Y 3 = l- Yi + 2X3T3 2 , ( 3 . 12 ) 

which is often called the autonomous limit. Since a generic integral curve given by Pih(,d 6 ) = 
A|T 3 2 — A 3 y 3 2 + A3 = const is an elliptic curve, a general solution can be expressed by Weier- 
strass’s elliptic functions. Then, the system (3.3) with small £3 can be studied by a perturbation 
method. 

Let us calculate the action of the extended affine Weyl group W ((2Ai)(P) restricted on the 
set {£3 = 0}, which leaves the autonomous limit (3.12) invariant. 

Proposition 3.4. The birational transformation group IT((2Ai)(P) is extended to the bira- 
tional transformation group acting on CP 3 (0,1, 2,1). The transformation group which leaves 
the autonomous limit (3.12) invariant is given by Z 2 k Aut ((2Ai)lP) = Z 2 k (th , 7 r2, 04 ) ■ 

Proof. The first part of Proposition is verified by a straightforward calculation. To show the 
second part, we should write down the actions in the (A3, Y 3 , £3)-chart. For example, the action 
of si in the (A3,13, £3)-chart is given by 

, £3, o, [5 J . 


(A3, Y3, £3, a, f 3 ) i-)- 


A 3 ,A 3 + a^ 
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On the set {£3 = 0}, it is reduced to 

{X 3 ,Y 3 ,a,/3)^(X 3 ,Y 3 ,-a,f3). 

Since the autonomous limit (3.12) is independent of the parameters a, (3, the action of si 
to (3.12) is trivial. Similarly, the actions of so, s' 0 and to (3.12) are reduced to the trivial one. 
On the other hand, it is easy to confirm that the restriction of the actions of 7ri, 7 t 2 , o~i are not 
trivial, which are explicitly given by 

7 T! : (X 3 , Y 3 ) ^ (X 3 Y 2 - Y 2 + 1, -I/F3), 
vr 2 : (X 3 ,Y 3 )^ (-X 3 Y£,1/Y 3 ), 

<?!■■ (X 3 ,Y 3 ) ^ (l-X; h V-lY 3 ). 


Furthermore, (3.12) is invariant under the Z 2 action Y 3 1— > —Y 3 due to the orbifold structure of 
CP 3 (0,1, 2,1). Since r = deg(z) = 2, (X 3 ,Y 3 , e 3 ) are coordinates on the lift C 3 of the quotient 
C 3 /Z 2 , where Z 2 action is defined by (X 3 ,Y 3 , £3) i->- (X3, — 13, — £3). This is reduced to the 
action Y 3 i->- — Y 3 on the set {£3 = 0}. ■ 

A similar result also holds for the other Painleve equations except for (Pyi) and summarized 
in Table 4. If r = deg(z) = 1, the symmetry group of the autonomous limit is given by the 
Dynkin automorphism group G = Aut ( R l 1 ') because the action of W (iT 1 ') on the set {£3 = 0} 
is reduced to the trivial action as above. If r > 1, the autonomous limit is further invariant 
under the Z r action arising from the orbifold structure of C P 3 (p, q, r, s). The autonomous limit 
on the set {£3 = 0} is not defined for (Pyi) in this way because r = 0. 



weight 

Rj 

symmetry 

Pi 

(3, 2,4, 5) 

X 2 - 4r 3 - 2 Y 

Z4 

Pii 

(2,1,2, 3) 

X 2 - Y 4 - Y 2 

Z 2 k Aut (A^) 

Piv 

(1,1,1,2) 

X 2 Y - XY 2 - 2XY 

Aut (A^ 1} ) 

P III(.D8) 

(-1,2,4,1) 

y‘2\x‘2 ly I 

AT 2 I 2Y 

Z4 K Z 2 

P III(D 7 ) 

(-1,2,3,1) 

X 2 Y 2 + X + Y 

Z3 k Aut (A^) 

PIII(Z3 6 ) 

(0,1, 2,1) 

X 2 Y 2 - XY 2 + A 

Z 2 k Aut ((2Ai)Pi) 

Pv 

(1,0,1,1) 

A 2 Y 2 - X 2 Y + AT 2 - AY 

Aut ( A 3 ] ) 


Table 4. Hamiltonian functions of the autonomous limit defined on the set {£3 = 0}. (A 3 , Y 3 ) is denoted 
by (A, Y) for simplicity. 

In the rest of this paper, the other Painleve equations (Pih(_d 7 ))! (Pm(D g ))> (Py) and (Pyi) 
are studied with the aid of the weighted projective spaces and dynamical systems theory (see [3] 
for (Pi), (Pii) and (Piy)). Since the strategy is completely the same as that for (Pih(d 6 )), we 
only show important steps and formulae. 


4 The third Painleve equation of type D 7 

The space CP 3 (—1,2,3,1) for (Phi(d 7 )) is covered by four inhomogeneous coordinates (Yi,Zi,£i), 
(A 2 ,Z 2 ,£ 2 ), (A 3 ,y 3 ,£ 3 ) and ( x,y,z ) related as 

x = £l = X 2 £ 2 = A 3 £ 3 , y = Yi£j “ 2 = £2 2 = Y 3 e^ 2 , z = ZiAf 3 = Z 2 e^ i = £j 3 . 
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We give the third Painleve equation of type D^ on the local chart ( x,y,z ). On the other 
local charts, (Pmm 7 )) is expressed as rational differential equations. By rewriting them as 
3-dimensional autonomous vector fields, we obtain 

Yi = 2Y 1 + 2Y? - Zi + ae\Y\, 

Z\ = 2>Z\ + 6 F 1 Z 1 — £iZ\ + 3asiZi, (4.1) 

£i = £i(l + 2Y\ + a£i), 

X.2 = 2 + 2X| — X 2 Z 2 + a£ 2 X 2 , 

Z 2 = &X 2 Z 2 + 3Z| + 3a:£2Z2 — 262^2, 

£2 = £2(2^2 + Z 2 + Q!£2), 

and 

X 3 - al 3 j , 

Y» + ah 3 j , (4.2) 



We have the decomposition 

CP 3 (—1, 2,3,1) = CP 2 (—1, 2,3) U C 3 (disjoint), 

where C 3 = {(x,y,z)} and CP 2 (—1,2, 3) = {ei = 0} U {£2 = 0} U {e 3 = 0} in coordinates. 
This implies that the set CP 2 (—1, 2, 3) is attached at infinity of the phase space C 3 . Thus, the 
asymptotic behavior of solutions can be studied by the limit £1 —> 0 or —> 0 or e 3 —> 0 . 

The autonomous limit is a Hamiltonian system obtained by putting e 3 = 0 for equation (4.2). 
On the set £ 3 = 0, the action of the extended affine Weyl group W(A^) ~ Aut(H^ 1 ' ) ) k W(A^) 
is reduced to the action of Aut(H[ ) ~ Z 2 defined by (A 3 ,y 3 ) i-)- {Y^,X^). Further, the 
autonomous limit is invariant under the Z 3 action (r = 3) induced by the orbifold structure, see 
Table 4. 

The vector field (4.1) has fixed points on the infinity set CP 2 (—1, 2, 3) given by (Yi, Zi, £ 1 ) = 
(—1,0,0) and (—1/2, —1/2,0). The Jacobi matrices of (4.1) at the fixed points are 

/ 2 1 cA /0 1 a/2 \ 

- 0 3 0 , - 3 0 (3a- l)/2 , 

\0 0 1 / \0 0 0 ) 


Xz = -1 - 2 x£Y 3 + £3 Q 
Y 3 = 1 + 2A 3 y 3 2 + £3 


respectively. The latter fixed point having a zero eigenvalue corresponds to the irregular singular 
point, and the former one corresponds to a movable singularity associated with the regular 
Laurent series solution (i) given in Section 2.2; The Laurent series (i) converges to the point 
(Yi, Zi, £ 1 ) = (—1, 0, 0) as z —> zq. The exceptional Laurent series solution (ii) does not converge 
to some point on CP 3 (—1, 2 , 3,1) as z —> zq. 

Remark 4.1. These fixed points are also included in the chart (X 2 >-^ 2 >£ 2 )- However, it is 
better to use the first chart (Y, Z[, £\) because the second chart has to be divided by the Z 2 
action due to the orbifold structure. 


To treat the Laurent series (ii), we use the Backlund transformation a defined by 
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We consider another space CP 3 (—1, 2, 3,1) with inhomogeneous coordinates denoted by (x, y, z) 
etc. We glue two copies of CP 3 (—1,2,3,1) by the transformation a. Then, it is easy to show 
that the exceptional Laurent series (ii) is converted to the regular series (i) in the (x, y, z )- chart, 
and it approaches to the fixed point (Yi, Zi, ej.) = (— 1 , 0 , 0 ) as z — > zq. 

To construct the space of initial conditions of (Pmm 7 )), we perform the weighted blow-ups 
at the points (Lj, Z\ , E\ ) = (—1,0,0) and (Lj, Z\ ,£\) = (—1,0,0). 

(i) blow-up at (Yi,Zi,ei) = (— 1 , 0 , 0 ). 

For the vector field (4.1), put Y\ =Y\ + 1 and 

u = Y\ — Z\ + aei, v = Z \, w = £ 2 - 


Then, the linear part is diagonalized around (—1,0,0) and we obtain 

ii = 2u — 2 u 2 + 2 uv + 4 v 2 — vw + a{uw — 2 to ), 
v = 3 v + v(w — 6u — 6u + 3 aw), 
w = w — w(2u + 2v — aw). 

The origin ( u , v, w ) = (0, 0, 0) is a singularity of the foliation of integral curves. To resolve it, 
we introduce the weighted blow-up defined by 

U = u\ = V 2 U2 = W3U3, V = u\v 1 = V'2 = W3V3, W = U\W\ = V2W2 = W3. 

The weight (2, 3,1) is taken from the eigenvalues of the Jacobi matrix at the singularity. The 
relation between the original coordinates ( x, y, z) and the new coordinates ( 113 , V 3 , W 3 ) is given by 

x = w 3 , y = u 3 + V3W3 - a/w 3 - l/w%, z = v 3 , 

U 3 = y — xz + a/x + 1/x 2 , W 3 = x, V 3 = z. ( 4 . 4 ) 


In the new coordinates, (Pih(d 7 )) is transformed into the Hamiltonian system 

dw3 dH\ du3 dHi 

dz du 3 ’ dz dw 3 ’ 

with the polynomial Hamiltonian function 

q ii 1 o n AO / O \ 

zH 1 = — u + u w — -w z + 2 uw z + w z — a[uw + w z) 


(4.5) 


(here, the subscript is omitted for simplicity). Furthermore, we can verify the symplectic rela¬ 
tion (3.7). 

(ii) blow-up at (Yi,Zi,ei) = (— 1 , 0 , 0 ). 

This singularity is resolved by the same way as above by the weighted blow-up with the weight 
(2,3,1). The result is easily obtained by the Backhand transformation a as follows. Define the 
new coordinates (uq,vq,wq) by 


uq = y — xz + a/x + 1 /x 2 , wq = x, vq = z. 


Substituting (4.3) yields 

u 6 = xz + xy - (1 - a)z/y + z 2 /y 2 , w 6 = -y/z, v 6 = -z. (4.6) 

By this coordinates change, (Piii(d 7 )) is transformed into the Hamiltonian system of (u6,we), 
whose Hamiltonian function has the same form as (4.5), for which a is replaced by a = 1 — a. 

In this manner, all singularities are resolved and the space of initial conditions of (Piii(d 7 )) 
is given by C 2 x ^ U C 2 U3 U C 2 ue w ^ glued by the symplectic transformations (4.4), (4.6). 
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5 The third Painleve equation of type D 8 

The orbifold CP 3 (—1, 2,4,1) for (Pih(d 8 )) is covered by four inhomogeneous coordinates related 
as 


x = £ 1 = x 2 £2 = X 3 £ 3 , y = Y 1 e 1 2 = £ 2 2 = Y 3 e 3 2 , Z = Zi£ 1 4 = Z 2 £ 2 4 = £g 4 . 


We give the third Painleve equation of type D$ on the local chart ( x,y,z ). On the other 
local charts, (Pmm 8 )) is expressed as rational differential equations. By rewriting them as 
3-dimensional autonomous vector fields, we obtain 


Yi = V? - 2F, 4 - YiZi, 

Zi = TY^Z X - 2 Z\ - 8YiZ\ + £l Y?Z U (5.1) 

X 3 = -2 X 3 Y 3 + - - ^2 + 4 £ 3^3, 

Y 3 = 2 X 3 Y 2 - y 3 Y 3 , (5.2) 

^ 2 

£ 3 = -4^3- 


We will not use the vector field written in (X 2 , Z 2 , £ 2 )-chart. 

The autonomous limit is a Hamiltonian system obtained by putting £3 = 0 for equation (5.2). 
It is known that (Pih(d 8 )) is invariant under the transformation n defined by 


(x,y,z) = ir{x,y,z) 


( xy 2 y z \ 


(5.3) 


On the set £3 = 0, this action is reduced to the Z 2 action given by (X 3 ,Y 3 ) i-)- (— X 3 Y 3 , 1 /I 3 ). 
Further, the autonomous limit is invariant under the Z 4 action (r = 4) induced by the orbifold 
structure, see Table 4. 

The vector field (5.1) has the fixed point (Yi, Zi, £ 1 ) = (1/2, 0, 0). The Jacobi matrix of (5.1) 
at the fixed point is 


1 

8 



The regular Laurent series solution (i) given in Section 2.2 converges to this point as z —>• zq, while 
the exceptional Laurent series solution (ii) does not converge to some point on CP 3 (—1, 2,4,1). 
To treat the Laurent series (ii), consider another space CP 3 (—1, 2,4,1) with inhomogeneous 
coordinates denoted by (ay, 2) etc. We glue two copies of CP 3 (—1, 2,4,1) by the transforma¬ 
tion (5.3). Then, it is easy to show that the exceptional Laurent series (ii) is converted to the reg¬ 
ular series (i) in the (x, y, i/- chart, and it approaches to the fixed point (Yi, Z\, £ 1 ) = (1/2, 0, 0) 
as z —> zo. 

To construct the space of initial conditions of (Pinm 8 )), we perform the weighted blow-ups 
at the points (Yi, Zi,£i) = (1/2,0,0) and (Yi, Z \, E \) = (1/2,0,0). 

(i) blow-up at (Yi,Zi,£i) = (1/2,0,0). 

For the vector field (5.1), put Yi = Yi — 1/2 and 


u = Yi — 2Zi, 


v = Zi, 


w = £ 2 ■ 
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Then, the linear part is diagonalized around (—1,0,0). To resolve the singularity ( u,v,w ) = 
(0,0,0) of the resultant equation, we introduce the weighted blow-up defined by 

U = u\ = V%U2 = W3U3, v = u\v \ = V2 = W3V3, W = U\W\ = V2W2 = W3. 

The relation between the original coordinates ( x,y,z ) and the new coordinates ( 113 ^ 3 , 103 ) is 
given by 


x = W3, 

2 1 

u 3 = y- 2 x z- 


y = U 3 + 2 v 3 wj + —2 , Z = V3, 

2 Wn 


W3 = X, 


v 3 = z. 


In the new coordinates, (Piii(_d 8 )) is transformed into the Hamiltonian system with the Hamil¬ 
tonian function 


zH\ = ^ + u 2 w 2 + w 2 z - | w 3 z + 4 uw A z + 4 w 6 z 2 


2 

W Z 


(here, the subscript is omitted for simplicity). Furthermore, we can verify the symplectic rela¬ 
tion (3.7). 

The blow-up at (Ij, Z 1 , ei) = (—1, 0, 0) is done in the same way and the Hamiltonian function 
is easily obtained by applying the transformation (5.3) to the above H\. In this manner, we can 
obtain the space of initial conditions. 


6 The fifth Painleve equation 

The orbifold CP 3 (1,0,1,1) for (Py) is covered by three inhomogeneous coordinates (lj,Zi,£i), 
(A 3 , 13 , 63 ) and ( x,y,z ) related as 

x = e^ 1 =X 3 £^ 1 , y = Y\ = Y3, z = Z 1 £^ 1 =£^ 1 . 

The second chart does not appear because q = 0. We give the fifth Painleve equation on the local 
chart ( x,y,z ). On the other local charts, (Py) is expressed as rational differential equations. 
By rewriting them as 3-dimensional autonomous vector fields, we obtain 

Yi = ~ 2 Yj + 21 2 — Y\Z\ + Y 2 Z\ + oq£ 1 — (oq + a;3)li£i, 

Z\ = —Z\ + 2I1Z1 — z 2 + £\Z\ + 21 1 Z~Y — (oq + a 3)ZiEi + Ol2£\Zi, 

£1 = £i(—1 + 21 1 — Z\ + 2A\Zi — {011 T Q^Gq T a^'Ziei); 

-^3 = -^3 + X 2 — 2X3Y3 — 2X 2 Y3 — Ol2£3 + («1 + 013 — l)£q A3, 

Yi = —13 + Y 2 — 2X3I3 + 2V3Y3 2 + ai£3 — (aq + aq)£3l3, 

£3 = —£ 3 . 

The autonomous limit is a Hamiltonian system obtained by putting £3 = 0 for equation (6.2). On 
the set {£3 = 0}, the action of the extended affine Weyl group W(A 3 ^) ~ Aut^^) k W{A^) 
is reduced to the action of Ard/A^) generated by (A 3 , 13 ) t-q (I 3 — 1 ,—A 3 ) and (A 3 ,l 3 ) eq 

(A 3 ,l -I3). 

The vector field (6.1) has fixed points on the infinity set CP 2 (1,0,1) given by 


( 6 . 1 ) 

( 6 . 2 ) 


{Yi,Zi,ei) = (0,0,0), (1,0,0), (0, -1,0), (1/2,-2,0), (1,-1,0). 
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The Jacobi matrices of (6.1) at these fixed points are 

/ 2 0 -aA / 2 0 -a 3 \ /-1 0 aA /0 -1/4 *\ 

- 01 0 , 01 0 , 0 1 * ) , [ 4 0 * 

\0 0 1 / \0 0 1 / \ 0 0 0 / \0 0 0 / 


0 

-1 

0 



respectively, where * denotes certain long numerical expressions. The latter three fixed points 
having zero eigenvalues correspond to the irregular singular point. The center-(un)stable mani¬ 
folds at these points can be exactly calculated for certain specific values of parameters, which give 
the Riccati equations of confluent hypergeometric type (Table 3). The former two points corre¬ 
spond to movable singularities associated with the regular Laurent series solutions (ii) and (i) 
given in Section 2.2. The Laurent series (ii) and (i) converge to the points (Yi, Z i, £i) = (0, 0, 0) 
and (1,0,0) as z —> zq, respectively. The exceptional Laurent series solutions (iii) and (iv) do 
not converge to some point on CP 3 (1, 0,1,1) as z — > zq. 

To treat the Laurent series (iii) and (iv), we use the Backlund transformation it defined by 

{x,y,z,ai,a 2 ,a 3 ) = n(x,y,z,ai,a 2 ,a 3 ) = [z{y - 1 ), z, a 2 , a 3 ,1 - aq - a 2 - a 3 ) . 

We consider another space CP 3 (1,0,1,1) with inhomogeneous coordinates denoted by (x, y, i) 
etc. We glue two copies of Ci- > 3 (1, 0,1,1) by the transformation n. Then, the exceptional 
Laurent series (iii) and (iv) are converted to the regular series (i) and (ii), respectively, in the 
(x, y, i)-chart. They converge to the fixed points (Y\,Z\,e\) = (1,0,0) and (0,0,0) as z —>• zq, 
respectively. 

To construct the space of initial conditions of (Py), we perform the weighted blow-ups at 
these four points. 

(i) blow-up at (Yi,Zi,ei) = (0,0,0). 

For the vector field (6.1), put u = Y\ — a\e i, v = Z\, w = sq. Then, the linear part is 
diagonalized. For the system of ( u,v,w ), we introduce the weighted blow-up defined by 

U = u\ = v\u 2 = Vj\u 3 , V = U\V\ = V 2 = W 3 V 3 , W = U\W\ = v 2 w 2 = w 3 . 

The weight (2,1,1) is taken from the eigenvalues of the Jacobi matrix at the singularity. The 
relation between the original coordinates (x, y, z) and the new coordinates (u 3 , v 3 , w 3 ) is given by 

x = 1/103, y = u 3 wj + aiw 3 , z = v 3 , 

u 3 = x 2 y — a\x. w 3 = 1/x, v 3 = z. 


In the new coordinates, (Py) is transformed into the Hamiltonian system, whose Hamiltonian 
function is 


zH\ = u — u 2 w 2 — u 2 w 3 z + (z — a\ + a 3 )uw — (2aq + a 2 )uw 2 z — a\{a\ + a 2 )wz 

(here, the subscript is omitted for simplicity). Furthermore, we can verify the symplectic rela¬ 
tion (3.7). 

(ii) blow-up at (Yi,Zi,ei) = ( 1 , 0 , 0 ). 

For the vector field (6.1), put Y\ = Y\ — 1 and u = Y\ — a 3 e 1 , v = Z\, w = s 1 . Then, the 
linear part is diagonalized around (1, 0, 0). For the system of (u, v, w), we introduce the weighted 
blow-up defined by 

U = u\ = V 2 Us = WqUq, V = U4X4 = U5 = WQVe, W = U4W4 = V5W5 = IUq. 

The relation between the original coordinates (x, y. z) and the coordinates (uq, vg,we) is given by 
x = 1 /w 6 , y = uqwI + a 3 WQ + 1, 


z = x 6 , 
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Uq = X 2 y — X 2 — CC3X, Wq = 1 /x, Vq = z. 

In the new coordinates, (Py) is transformed into the Hamiltonian system, whose Hamiltonian 
function is 

zH\ = —u — u 2 w 2 — u 2 w 3 z + (—z + a\ — a^)uw — ( 2 a 3 + a.2)uw 2 z — 03(0:2 + a^)wz 

(here, the subscript is omitted for simplicity). Again, we can verify the symplectic relation (3.7). 

The singularities (Yl,Zi,£i) = (1,0,0) and (0,0,0) are resolved by the same way as above 
by the weighted blow-up with the weight (2,1,1). The result is easily obtained by the Backlund 
transformation it as in the previous sections. In this manner, all singularities are resolved and 
it turns out that the space of initial conditions of (Py) is given by five copies of C 2 glued by the 
symplectic transformations. 

7 The sixth Painleve equation 

The orbifold CP 3 (1, 0,0,1) for (Pyi) is covered by two inhomogeneous coordinates (Yi,Zi,£i) 
and ( x , y , z) related as 

x = ej" 1 , y = Yi, z = Zi. 

We give the sixth Painleve equation on the local chart (x,y,z). On the other local chart 
(Yi,Zi,ei), (Pyi) is expressed as a rational differential equation. As before, we rewrite it as 
a 3-dimensional autonomous polynomial vector field, whose expression is too long and omitted 
here. This vector field has fixed points 

(Yi,Zi,ei) = (0,c,0), (1,c,0), (c,c,0), 

where c £ C is an arbitrary constant. The Jacobi matrices at these fixed points are 

/2 0 — a 4 \ /2 0 — a 3 \ /2 -2 1 - a 0 \ 

c 0 0 c-1 , (1-c) 0 0 -c , c(c-l) 0 0 1 , 

\0 0 1 / \0 0 1 / \0 0 1 / 

respectively. These fixed points correspond to movable singularities associated with the regular 
Laurent series solutions (iii), (ii) and (i) given in Section 2.2. The Laurent series (iii), (ii) and (i) 
converge to these points as z —> zo, respectively. The exceptional Laurent series solutions (iv) 
and (v) do not converge to some point on CP 3 (1,0, 0,1) as z —> zo. 

To treat the Laurent series (iv) and (v), we use the Backlund transformation tt\ defined by 

(■E U1 ~i ^0; Q- 1 1 0-2 ; ^3; ^4) ^"l ( 4' ■ 2/, Zj Oi 0 , Oil, CK2 ; 0:3; ^4) 

= ( ~~{xy + a 2 ), ~,z,a 3 , 04,0:2, ao,ai 

V ^ y 

We consider another space CP 3 (1,0,0,1) with inhomogeneous coordinates denoted by (x, y, z) 
etc. We glue two copies of CP 3 (1,0,0,1) by the transformation tt\. Then, the exceptional 
Laurent series (iv) becomes the regular series (iii) in the (x, y, 5)-chart. It converges to the fixed 
point (Yi,Zi,£i) = (0, c, 0) as z —> zq. On the other hand, the exceptional Laurent series (v) is 
expressed as x(z) o(i), y(z) 0(T ) in the (x, y, z)-chart. Since the solution is holomorphic 
at z = zo, we need not resolve the singularity. 

To construct the space of initial conditions of (Pyi), we perform the weighted blow-ups at 
four points. 
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(i) blow-up at (Yi,Zi,ei) = (0,c, 0). 

For the vector field written in (Yi, Z\, £i)-chart, put u = Y\ — a.±£\, v = Z\, w = £1. Then, 
we introduce the weighted blow-up with the weight (2, 0,1) defined by 

U = lt\ = W2U2, V = V\ = V 2 , W = U\W\ = VJ2- 

This is a blow-up along the line {v = const}. The relation between the original coordinates 
(, x,y,z ) and the new coordinates (u2,V2,W2) is given by 

x = l/w2, y = U2V0I + a 2 W2 } z = v 2 , 

U2 = x 2 y — a.2X, W2 = 1/a;, V2 = z. 

In the new coordinates, (Pyi) is transformed into the Hamiltonian system, whose Hamiltonian 
function is 

z(z — l)Hi = u 3 w 4 + uz — u 2 w 2 ( 1 + z) + (1 — «o — «3 + 2a.4)u 2 w 3 

+ (ao — «4 + za-i — za.4 — 1 )uw + 04(0102 + of + (1 — 00 — 03)04)?/; 

T (01O2 + O2 2 a 4 — 200O4 — 203O4 T Oj) // /r 

(here, the subscript is omitted for simplicity). Furthermore, we can verify the symplectic rela¬ 
tion ( 3 . 7 ). 

The other singularities (Yi,Zi,£i) = ( 1 , c, 0 ), (c, c, 0 ) and (Yi,Zi,£i) = ( 0 ,c, 0 ) are resolved 
by the same way as above by the weighted blow-up with the weight ( 2 , 0 , 1 ). In this manner, all 
singularities are resolved and the space of initial conditions of (Pyi) is obtained by six copies 
of C 2 glued by the symplectic transformations; i.e., ( x,y ), (x,y) and four charts arising from 
blow-ups at four points. We need (x,y )-chart for the exceptional Laurent series (v). 
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